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DERIVATION  OF  TRANSFER  FUNCTIONS  FOR  A  FLUID- 
LOADED,  MULTIPL-LAYER  THICK  PLATE  SYSTEM 


1.  INTRODUCTION 


Transfer  functions  are  mathematical  expressions  that  relate  two  physical  quantities  of  a 
system.  When  applied  to  mechanical  systems,  transfer  functions  model  energy  propagation 
through  the  medium  of  a  structure.  This  medium  can  be  simple  or  complex,  depending  on  the 
type  and  composition  of  the  structure.  Typically,  partial  differential  equations  describe  some 
field  quantity  on  the  medium.  A  set  of  assumptions  is  made,  and  these  partial  differential 
equations  reduce  to  ordinary  differential  equations  and  then  to  algebraic  equations.  The  end 
result  is  an  expression  that  is  the  quotient  of  two  of  the  field  quantities,  frequently  a  system  value 
at  some  location  in  the  field  divided  by  an  energy  input  into  the  field.  These  transfer  functions 
are  useful  in  understanding  system  behavior,  and  they  sometimes  can  be  verified  by 
measurements  at  different  locations  in  the  structure. 

This  report  derives  six  mathematical  transfer  functions  of  various  thick  plate  systems. 
First,  two-dimensional  equations  of  elasticity  are  derived  and  mathematical  expressions  for  the 
plate  displacements  and  spatial  derivatives  of  the  plate  displacements  are  calculated.  Second,  the 
equations  of  motion  of  a  single  thick  plate  with  a  normal  pressure  load  are  then  determined. 
Third,  the  single  plate  is  coupled  with  a  fluid  whose  dynamics  are  governed  by  the  wave 
equation.  This  fluid  contains  an  acoustic  pressure  that  acts  on  the  structure  and  changes  the 
dynamic  behavior  of  the  system.  The  response  of  this  system  is  calculated.  Fourth,  the 
equations  of  motion  of  two  connected  thick  plates  with  a  normal  pressure  load  are  determined. 
This  double  (thick)  plate  allows  for  two  different  types  of  plate  materials  to  be  attached  together. 
Fifth,  the  fluid  load  is  then  coupled  with  the  double  thick  plate  and  the  response  is  calculated. 
Sixth,  the  equations  of  motion  of  three  connected  thick  plates  with  a  normal  pressure  load  are 
determined.  The  triple  (thick)  plate  allows  for  three  different  types  of  plate  materials  to  be 
attached  together.  Seventh,  the  fluid  load  is  then  coupled  to  the  triple  thick  plate  and  the 
response  is  calculated.  In  all  cases,  the  thick  plate  equations  of  motion  are  verified  by  comparing 
their  transfer  functions  to  a  second  transfer  function  derived  using  the  thin  plate  theory  calculated 
with  a  very  small  plate  thickness.  This  allows  validation  of  the  results  using  a  second  method  of 
modeling  the  response. 


2.  EQUATIONS  OF  ELASTICITY 


The  transfer  functions  are  derived  by  assuming  that  each  thick  plate  is  governed  by  the 
equation  (reference  1) 


2  <7  u 

p57 u  +  (A  +  /t)VV  •  u  =  p — Y  ■ 

dt 


(1) 
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where  p  is  the  density;  A  and  p  are  the  Lame  constants;  t  is  time;  •  denotes  a  vector  dot  product; 
u  is  the  Cartesian  coordinate  displacement  vector  expressed  as 


ux(x>y>z,t) 

u=-|wy(^,y,z,OK 

uz(x,y,z,t)_ 


(2) 


with  subscript  x  denoting  the  direction  parallel  to  the  plates,  y  denoting  the  direction  into  the 
plates,  and  z  denoting  the  direction  perpendicular  to  the  plates.  The  modeled  geometry  and  the 
coordinate  system  of  the  plates  are  shown  in  figure  1.  The  symbol  V  is  the  gradient  vector 
differential  operator  written  in  three-dimensional  Cartesian  coordinates  (reference  2)  as 


V 


d  .  1  d  .  d  . 

frlr+7!)eld+!klz' 


(3) 


with  ix  denoting  the  unit  vector  in  the  ^-direction,  i  denoting  the  unit  vector  in  the  y-direction, 

and  iz  denoting  the  unit  vector  in  the  z-direction;  V2  is  the  three-dimensional  Laplace  operator 
operating  on  vector  u  as 


V2u  =  V2uxix  +  V2uyiy  +  V2u  X  , 


with  V2 


operating  on  scalar  u  as 


(4) 


Ux,y<z=V*^Ux,y,: 


d\y.z  d_\±z 

dx 2  dy 2  dz 2 


and  the  term  V  *u  is  called  the  divergence  and  is  equal  to 


V»u  = 


duy 

dy 


(5) 


(6) 


The  applied  loading  effects  of  the  pressure  in  the  absence  of  a  fluid  acting  on  the  surface 
of  the  plate  are  modeled  as  a  structural  load  in  sections  3, 5,  and  7.  The  applied  loading  effects 
of  the  pressure  in  the  presence  of  a  fluid  acting  on  the  surface  of  the  plate  are  modeled  as  an 

acoustic  load  in  sections  4,  6,  and  8.  Note  that  p.  A,  and  u  are  the  properties  of  each  specific 
plate. 


The  displacement  vector  u  is  written  as 
u  =  V0  +  Vxt/?  ? 
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where  <p  is  a  dilatational  scalar  potential,  x  denotes  a  vector  cross  product,  and  V  is  an 
equivoluminal  vector  potential  expressed  as 


¥x(x,y,z,t) 
\jf  =  <y/y(x,y,z,t) 
y/z(x,y,z,t) 


(8) 


The  problem  is  now  formulated  as  a  two-dimensional  response  (y  =  0  and  d(-)/dy  =  0)  problem. 
Expanding  equation  (7)  and  breaking  the  displacement  vector  into  its  individual  nonzero  terms 
yield 


_d$(,x,z,t)  dVy&zS) 

dx  dz 


(9) 


and 

d(j)(x,z,t)  dvy(x,z,t) 

u,(x,z,t)  = - - - +  — - 

dz  dx 


(10) 
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Equations  (9)  and  (10)  are  next  inserted  into  equation  (1),  which  results  in 


c2dV2(j){x,z,t ) 


<?2(ft(x,z,0 

dt 2 


and 


c2vV/x,z,r) 


dVyCW) 

dt 2 


(ID 


(12) 


on  all  material  layers.  Equation  (1 1)  is  the  dilatational  component  and  equation  (12)  is  the  shear 
component  of  the  displacement  field  (reference  3).  Correspondingly,  the  constants  cd  and  cs  are 

the  complex  dilatational  and  shear  wave  speeds,  respectively,  and  are  determined  by 


cd 


(13) 


and 


(14) 


The  relationship  of  the  Lame  constants  to  the  compressional  and  shear  moduli  is  shown  as 


and 


A  = 


Ev 

(1  +  i>)(l  -  2v)  ’ 


fi  -  G  = 


E 

2(l+u)’ 


(15) 

(16) 


where  E  is  the  complex  compressional  modulus  (N/m^),  G  is  the  complex  shear  modulus  (N/m^), 
and  v  is  the  Poisson's  ratio  of  the  material  (dimensionless). 

The  conditions  of  infinite  length,  two-dimensional  response  (y  =  0),  and  steady-state 
response  are  now  imposed,  allowing  the  scalar  and  vector  potential  to  be  written  as 

<j)(x,z,t)  =  <I>(z)exp(i^;cx)exp(ifyr)5 
and 

y/y(x,z,t)  =  xF(z)exp(iLcx)exp(iftW) , 

where  kx  is  the  wavenumber  of  excitation  in  the  x-direction  (rad/m),  co  is  the  frequency  of 

excitation  (rad/s),  and  i  is  the  square  root  of  -1.  Note  that  equations  (17)  and  (18)  are  valid  on 
every  layer. 

Inserting  equation  (17)  into  equation  (11)  yields 


(17) 

(18) 
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d2*(z)  ,  _,2 


dz 2 


+a  <J>(z)  =  0, 


where 


and 


a  =  Jkf-ki . 


kd  = 


co 


(19) 


(20) 


(21) 


Inserting  equation  (18)  into  equation  (12)  yields 

<f2*U)  ,  *2, 


dz 2 


■+j0^(z)  =  O, 


where 


and 


cr 


The  solution  to  equation  (19)  is 
4>(z)  =  A(^,ft))exp(ioe)  +  5(/:x,co)exp(-iaz) , 
and  the  solution  to  equation  (22)  is 

'P(z)  =  C(kx,(o)exp(i/3z)  +  £>(fcx,co)exp(-i/?z) , 


(22) 


(23) 


(24) 


(25) 


(26) 


where  A,  B,  C,  and  D  are  constants  that  are  determined  below.  The  displacements  can  now  be 
written  as  functions  of  the  unknown  constants.  They  are 

uz (x, z, t)  =  {  i a[A{kx , co) exp(i az)-B(kx,(0 ) exp(-i ccz)]+ 

ikx[C(kx,0))exp(i(3z)  +  D(kx,a))exip(-ifiz)]  }exp(i^xx)exp(icot) ,  (27) 


and 


ux(x,z,t )  =  {  i£x[A(fcx,co)exp(ic#)  +  i?(fc.t,co)exp(-icc)]- 

iP[C(kx,a))txp(iPz)-D(kx,co)exp(-iPz)]  Jexp(i^xx)exp(icor) .  (28) 
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The  solution  to  the  constants  is  determined  by  formulating  the  problem  based  on  the  number  of 
plates  and  the  presence  or  absence  of  a  fluid  load.  This  is  done  in  sections  3  -  8.  Four 
derivatives  of  equations  (27)  and  (28),  which  are  used  in  the  solution  of  the  constants,  are 

dur(x,z,t)  r  ,?r  n 

- ^ ~ kx  [AfA^colexpOa?)  +  B(kx,(o)exp(-iaz)]+ 

fax [C(kx,a))exp(ipz) - D(kx,(o)exp(-iPz)]  }exp(i^x)exp(i<yr) ,  (29) 

dur{x,z,t)  r  ,  r  i 

- ^ - =  I  - kxa[A(kx,(0)exp(iaz) - B(fcJC,(o)exp(-iaz)]+ 

P2[C(kx,Q))exip(iPz)  +  D(kx,a>)exTp(-i(3z)]  }exp(i£xx)exp(iatf),  (30) 


— ^  =  {  -«2[A(/:;c,£y)exp(iaz)  +  5(A:;(.,(w)exp(-iaz)]- 

fax  ic (kx , e*p(i/fe)  ~D(kx, (o) exp(-i/fe)]  }exp(ifcxx)exp(ia*) , 

— ~  i  kxalMkx , cd) exp(i az)  -  5(^,m)exp(-iaz)]- 

kx[C(kx,CQ)exp(iPz)  +  D(kx,(D)exp(-ifiz)]  }exp(iA:;(.x)exp(ift)r) 


(31) 


(32) 
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3.  SINGLE  PATE  WITH  NO  FLUID  LOAD 


The  first  transfer  function  derived  is  a  sigle  plate  with  no  fluid  load.  The  structural  load 
consists  of  an  applied  pressure  at  definite  wavenumber  and  frequency  at  location  z  =  b,  as  shown 
in  figure  2.  The  normal  and  tangential  stresses  in  the  system  at  the  boundaries  are 


rzz(x,b,t)  =  (X  +  2p) 


duz(x,b,t)  +^dux(x,b,t) 


dz 


dx 


=  -ps(xM > 


Ta(x,b,t)  =  n  j 


dux(x,b,t )  duz(x,b,t ) 
dz  dx 


=  0 


(33) 

(34) 


z.Jx,a,t)  =  (X  +  2p) 


du,(x,a,t )  dux(x,a,t ) 


* 


+  A- 


=  0, 


(35) 


and 


(^»^>0  —  P 


dux(x,a,t )  du  (x,a,t) 


dz 


dx 


=  0 


(36) 


where  ps(x,b,t )  is  the  structural  load  applied  at  z  =  b.  This  load  is  modeled  as  a  function  at 


definite  wavenumber  and  frequency  as 

ps(x,b,t )  =  exp(i^x)exp(ift») .  (37) 

Combining  equations  (29)  -  (37)  yields  the  four-by-four  linear  system  of  equations 

Ax  =  b,  (38) 

where  the  entries  of  equation  (38)  are 

Aj  j  =  (-  a2X  -  2a1  fi  -  Xkx  )exp(i  ab ) ,  (39) 

A12  =  (- cc2X-2a2p-Xkx)exp(-iab ) ,  (40) 

Aj  3  =  (-  2kx/3fi)exp(ifib) ,  (41) 

ai,4  =  (2kxPp)zxp(-iPb) ,  (42) 

A24  =  (-  2pkxa)exp(i  ab) ,  (43) 

A2  2  =  (2pkxa)exp(-iab) ,  (44) 
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4,3  =  (u/32  -^)exp(i^). 

(45) 

4,4  ={^P2~^ )exp(-i fib) , 

(46) 

4,i  =  (“  a2 &  -2a2  jx-  Afc2  )exp(iaa) , 

(47) 

4,2  =  (“«2^-  -  2  a2n~Xkl  )exp(-iaa) , 

(48) 

4,3  =  (“  2kxfi/j,)txp(ifia) , 

(49) 

4,4  =  (24^)exp(-i  fid) , 

(50) 

4i  =  (_  2/dkxa)exp(\aa) , 

(51) 

4,2  =  (2fJkxa)txp(-icca) , 

(52) 

4,3  =  ifiP2  ~  M2)ex p(ifia) , 

(53) 

4.4  =  (u£2  -  pLk\  )exp(-i  fia) , 

(54) 
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*u  =A(kx,a>),  (55) 

x2A  =  B(kx,co) ,  (56) 

x3l=C(kx,co),  (57) 

x41=D(kx,o)),  (58) 

hi  =-Ps(kx,a>),  (59) 

*2,1  -0,  (60) 

hi=0’  (6i) 

and 

*4j=0.  (62) 

Using  equations  (39)  -  (62),  the  solution  to  the  constants  A,  B,  C,  and  D  can  be  found  by 

x  =  A_1b.  (63) 


Additionally,  the  transfer  function  between  the  normal  displacement  at  location  z  and  the 
structural  load  can  be  written  as 


T(kx,(o)  =  =  A(kx,(o)iaexp(iccz)  -  B(kx,co)ia&xp(-iaz)  + 

rs  (Kx  ,(D) 

C(kx,a))ikx  exp(i/fe)  +  D{kx,(o)\kx  exp(-i/%) , 


or  the  normal  stress  and  the  structural  load 


T  (kx  ,0) )  =  =  Mkx  ,o>)(-  a2  A  -2a2  fi-  Xk 2  )exp(iaz)  + 

Ps  (kx ,  CO) 

B(kx,co)( -  a2 X  -  2a2 n  -Xkx  )exp(-iaz)  + 

C(kx ,  (Q)(-2kxpn)  exp(i/fe)  +  D(kx  ,0)){2kxfiix)txg{-\fiz) .  (65) 

The  transfer  function  in  equation  (64)  can  be  compared  to  another  transfer  function 
derived  using  the  thin  plate  theory  to  ensure  that  it  is  accurate.  The  thin  plate  equation  of  motion 
is 
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(66) 


-d4u  (x,t)  d2u.(x,t ) 

B~d? —  ph — P —  =  ~Ps  {x,t)  ’ 

where  B  is  the  bending  stress  expressed  as 


B  = 


Eh3 


12(1  -vl) 


(67) 


and  h  is  the  thickness  of  the  plate  (m)  and  corresponds  to  b-a  in  figure  2.  The  transfer  function 
between  transverse  displacement  and  applied  force  for  the  thin  plate  is 


T(k  ,(0)  = 


UJkx,co) 


-1 


Ps(kx,(o)  (Bkx  -  ph(02)' 
The  flexural  wavenumber  of  this  system  is  at 


phor 

B 


-il/4 


(68) 


(69) 


where  kj  has  units  of  rad/m.  At  this  wavenumber,  the  response  of  the  system  will  be  extremely 
large  compared  to  its  response  at  other  wavenumbers. 

Figure  3  is  a  plot  of  displacement  in  the  ^-direction  (normal)  divided  by  normal  pressure 
versus  wavenumber  in  the  ^-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  =  b.  The 
solid  line  represents  the  thick  plate  theory  (equation  (64))  and  the  x’s  represent  the  thin  plate 
theory  (equation  (68)).  The  flexural  wavenumber  calculated  from  equation  (69)  is  denoted  on  the 
plot.  In  this  example,  the  thickness  of  the  plate  is  very  small  (h  =  0.01  m)  so  that  a  direct 
comparison  between  the  thick  and  thin  plate  theories  can  be  made.  Other  parameters  used  to 
formulate  this  model  are:  bottom  of  the  plate  (a)  is  -0.01  m,  top  of  the  plate  (b)  is  0  m.  Young’s 

modulus  ( E )  is  13.2e9  (l+0.03i)  N/m2,  Poisson’s  ratio  (u)  is  0.30,  and  density  (p)  is  1938 
kg/m  .  Only  positive  values  of  wavenumber  are  shown,  as  the  function  is  symmetric  about  k  = 

0.  Note  that  there  is  almost  complete  agreement  between  the  thin  and  thick  plate  theories  for  this 
specific  example.  Figure  4  is  a  plot  of  displacement  in  the  z-direction  (normal)  divided  by 
normal  pressure  versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  and  a 
plate  thickness  of  h  =  0.1  m  evaluated  at  z  =  b.  The  two  models  begin  to  diverge  as  the 
wavenumber  increases.  This  is  due  to  the  inclusion  of  the  rotary  inertia  terms  in  the  thick  plate 
model  that  are  not  included  in  the  thin  plate  model.  Figure  5  is  a  plot  of  the  wavenumber- 
frequency  (kco)  surface  of  the  system  using  the  thick  plate  equations  of  motion  (equation  (64)) 
and  the  parameters  that  were  used  to  formulate  figure  4.  The  color  scale  to  the  right  of  the  plot  is 
in  decibels. 
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Phase  Angle  (degrees)  Magnitude  (dB) 


Figure  3.  Transfer  Function  of  a  Single  Plate  with  No  Fluid  Load 
( Equations  (64)  and  (68))  at  500  Hz  with  h  =  0.01  m 
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Figure  4.  Transfer  Function  of  a  Single  Plate  with  No  Fluid  Load 
(Equations  (64)  and  (68))  at  500  Hz  with  h-0.1m 


Figure  5.  Transfer  Surface  of  a  Single  Plate  with  No  Fluid  Load 
(Equation  (64))  with  h  =  0.1  m 


13/(14  blank) 


4.  SINGLE  PLATE  WITH  A  FLUID  LOAD 


The  second  transfer  function  derived  is  a  single  plate  with  a  fluid  load  on  one  side.  The 
applied  load  is  modeled  as  an  incident  pressure  wave  in  the  fluid  at  definite  wavenumber  and 
frequency  at  location  z  =  b,  as  shown  in  figure  6.  The  normal  and  tangential  stresses  in  the 
system  at  the  boundaries  are 


oz  ox 


Tzx(x,b,t)  =  fi 


dux(x,b,t)  du.(x,b,t ) 


dz 


dx 


=  0 


(70) 

(71) 


dz 


dx 


and 


rzx(x,a,t)  =  n 


dux(x,a,t )  du  (x,a,t) 


dz 


+  - 


dx 


=  0. 


(72) 


(73) 


where  pa(x,b,t)  represents  the  (acoustic)  pressure  of  the  fluid  load  on  the  plate  and  includes  the 
applied  load. 

The  acoustic  pressure  in  the  outer  fluid  is  governed  by  the  wave  equation  and  is  written  in 
Cartesian  coordinates  (reference  4)  as 

<?2j>aO,Z,t)  <?2R,(.t,Z,<)  1  <?2/>„(jt,Z,<) 

A2  A2  c)  dt1  (  ' 

where  pa(x,z,t)  is  the  pressure  ( N/m2),  z  is  the  spatial  location  (m)  normal  to  the  plate,  and  c  f 

is  the  compressional  wavespeed  of  the  fluid  (m/s).  The  acoustic  pressure  is  modeled  as  a 
function  at  definite  wavenumber  and  frequency  as 

Pa  (x,  z,t)  =  Pa  (z,kx,co)exp(ikxx)exp(ia)t) .  (75) 

Inserting  equation  (75)  into  (74)  and  solving  the  resulting  ordinary  differential  equation  yields 

Pa(z,kx,(o)  =  G(kx,o})exp(iyz)  +  PI(kx,(o)e\p(-iyz ),  (76) 

where  the  first  term  on  the  right-hand  side  represents  the  reradiated  pressure  field  and  the  second 
term  represents  the  applied  incident  pressure  field  (the  forcing  function)  acting  on  the  structure. 
In  equation  (76), 
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Fluid 


Figure  6.  Single  Plate  with  Fluid  Load 


r  = 


r  \2 

ft) 


-K 


(77) 


where  yis  purely  real  or  imaginary,  depending  on  the  sign  of  the  argument.  When  the  sign  of  the 
argument  is  positive,  the  analysis  is  in  the  acoustic  cone;  when  the  sign  of  the  argument  is 
negative,  the  analysis  is  in  the  nonacoustic  region.  The  interface  between  the  fluid  and  solid 
surface  of  the  plate  satisfies  the  linear  momentum  equation  (reference  5),  which  is 

d2uz(x,b,t)  dp  a  (x,b,t) 

Pf  ~2  -  — - ,  (78) 

where  Pj  is  the  density  of  the  fluid  (kg/m3 ). 

Combining  equations  (70)  -  (78)  yields  the  four-by-four  linear  system  of  equations 
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Ax  =  b , 


(79) 


where  the  entries  of  equation  (79)  are 


4,i  = 


2^ 


-a2A-2a2/t-/U:2  +  Pf03  a 


exp(iafr) , 


J 


4,2  ~ 


'  2,  ,_2..  1 1.2  Pr<»2a' 


-a X-2a p-Ak2 -■ 


v 


exp(-i  ab). 


4,3  “ 


-2  kxPn  + 


pf(o\^ 


exp(i fib) , 


4,4  - 


2kxfip  + 


Pf(02kx^ 


exp(-i  fib), 


\  J 

4.1  =  (~2pkxa)exp(iab), 

4.2  =  (2pkxa)exp(-iab) , 

4.3  =  (pP2~Pkx  )exp(ij3b) , 

4.4  =  (pP2  ~ pk2)exp(-ifib) , 

4.1  =  (-  ~  2a2 p  -  Ak2x  )exp(ia<2) , 

4.2  =  (-  vPh  ~  2 a2  p  - Ak2x  )exp(-iaa) , 

4.3  =  (-  2kxfip)exp(ifia) , 

4.4  =  (2kxfip)exp(-ifia) , 
A41=(-2pkxa)exp(icxa), 

4.2  =  (2pkxa)exp(-iaa) , 

4.3  =  (pP2  ~  Pkx  )exp(i/Ja) , 


(80) 

(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 
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(95) 


A4,4  =  (j*PZ  -/ifc2)exp(-i pa), 
xu=A(kx,0)),  (96) 

x2l=B(kx,(o),  (97) 

*3,i  =C(kx,co),  (98) 

xA1=D(kx,co),  (99) 

bll=-2PI(kx,(o)exp(-iyb),  (100) 

^2,1  =  0  .  (101) 

^3,1=0,  (102) 

and 

^4,i=0.  (103) 

Using  equations  (80)  -  (103),  the  solution  to  the  constants  A,  B,  C,  and  D  can  be  found  by 

x  =  A_1b.  (104) 

Additionally,  the  transfer  function  between  the  normal  displacement  at  location  z  and  the 
acoustic  load  can  be  written  as 

T(kx,Q) )  =  - - ~  =  A{kx,(o)'\atxp{}(xz)-B{kx,(o)\atx^{-iaz)  + 
P,{kx,w)exv(-\yb) 

C(kx,co)ikx  exp(i fk)  +  D(kx,(o)ikx  exp(-i/%) ;  qo5) 

or  the  normal  stress  and  the  structural  load 

T(K,a»  =  =  A(t„ffl)(-a!A  -  2o  v  -  Ai’)exp(i<K)  + 

P!(kx,(o)txp(-iyb) 

B( kx ,  to)(-  a 2 A  -  2a 2 pi  - Xk2x  )exp(-i az)  + 

C  ( kx ,  (0)(-2kxpn )  exp(i  /fe)  +  D{kx ,  c o){2kx  f$n)  exp(-i/fe).  ( 1 J6" 

The  transfer  function  in  equation  (105)  can  be  compared  to  another  transfer  function 
derived  using  the  thin  plate  theory  to  ensure  that  it  is  accurate.  This  transfer  function  is  derived 
by  equating  equation  (66)  to  equation  (76)  and  applying  equation  (78),  which  yields 
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(107) 


(108) 


t 

where  k  y  has  units  of  rad/m  and  £y  is  determined  using  equation  (69).  At  this  wavenumber,  the 
response  of  the  system  will  be  extremely  large  compared  to  its  response  at  other  wavenumbers. 

Figure  7  is  a  plot  of  displacement  in  the  z-direction  (normal)  divided  by  normal  pressure 
versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z-b.  The 
solid  line  represents  the  thick  plate  theory  with  a  fluid  load  (equation  (105))  and  the  x’s  represent 
the  thin  plate  theory  with  a  fluid  load  (equation  (107)).  The  fluid-loaded  flexural  wavenumber 
calculated  from  equation  (108)  is  denoted  on  the  plot.  In  this  example,  the  thickness  of  the  plate 
is  very  small  (h  =  0.01  m)  so  that  a  direct  comparison  between  the  thick  and  thin  plate  theories 
can  be  made.  Other  parameters  used  to  formulate  this  model  are:  bottom  of  the  plate  (a)  is  -0.01 

m,  top  of  the  plate  ( b )  is  0  m.  Young’s  modulus  (E)  is  13.2e9  (l+0.03i)  N/m2,  Poisson’s  ratio 
(u)  is  0.30,  density  of  the  plate  (p)  is  1938  kg/m3 ,  compressional  wavespeed  of  the  fluid  (cy)  is 

1500  m/s,  and  density  of  the  fluid  (py)  is  1025  kg/m3 .  Only  positive  values  of  wavenumber  are 

shown,  as  the  function  is  symmetric  about  k  =  0.  Note  that  there  is  almost  complete  agreement 
between  the  thin  and  thick  plate  theories  for  this  specific  example.  Figure  8  is  a  plot  of 
displacement  in  the  z-direction  (normal)  divided  by  normal  pressure  versus  wavenumber  in  the  x- 
direction  at  a  forcing  frequency  of  500  Hz  and  a  plate  thickness  of  h  =  0.1  m  evaluated  at  z  =  b. 
The  two  models  begin  to  diverge  as  wavenumber  increases.  This  is  due  to  the  inclusion  of  the 
rotary  inertia  terms  in  the  thick  plate  model  which  are  not  included  in  the  thin  plate  model. 

Figure  9  is  a  plot  of  the  wavenumber-frequency  (kco)  surface  of  the  system  using  the  thick  plate 
equations  of  motion  (equation  (105))  and  the  parameters  that  were  used  to  formulate  figure  7. 

The  color  scale  to  the  right  of  the  plot  is  in  decibels. 
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Figure  7.  Transfer  Function  of  a  Single  Plate  with  Fluid  Load 
( Equations  (105)  and  (107))  at  500  Hz  with  h  =  0.01  m 
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Figure  8.  Transfer  Function  of  a  Single  Plate  with  Fluid  Load 
( Equations  (105)  and  (107))  at  500  Hz  with  h  -  0.1  m 
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Figure  9.  Transfer  Surface  of  a  Single  Plate  with  Fluid  Load 
(Equation  (105))  with  h  =  0.1  m 


5.  DOUBLE  PLATE  WITH  NO  FLUID  LOAD 


The  third  transfer  function  derived  is  a  double  plate  with  no  fluid  load.  The  structural 
load  consists  of  an  applied  pressure  at  definite  wavenumber  and  frequency  at  location  z  =  c,  as 
shown  in  figure  10.  The  normal  and  tangential  stresses  in  the  system  at  the  boundary  z  =  c  are 


/  *  *  .  du„(x,c,t )  ,  dur(x,c,t ) 

^ zz  (X’C’0  —  (^2  +  2/^2  )  i~T'  T  -  —  Ps  (x»CtO> 


dz 


dx 


(109) 


and 


^ zx  (x>  ^’0  ~ 


dux(x,c,t)  +  duz(x,c,t) 


dz 


dx 


=  0. 


(110) 


where  the  subscript  2  denotes  plate  2.  The  interface  between  plates  2  and  1  requires  four 
equations.  The  first  two  are  displacement  constraints,  which  are 

■  an) 

and 

u^bAv^2  =  ux^M\^  .  (112) 


The  second  two  are  stress  constraints,  which  are 

Tzz  (*>  b,  0|plate  2  =  Tzz  (x’  b,  0|platel  > 

and 

T^^fc’r)|plate2=T-^i,’f)|platel  * 


(113) 

(114) 


Finally,  the  normal  and  tangential  stresses  in  the  system  at  the  boundary  z  =  a  are 


dz 


dx 


and 


,  .  r  dux(x,a,t)  duz(x,a,t)~\ 

Tzx(x,a,t)  =  m\ - - - +  - 1  =  °  > 


(115) 


(116) 


where  the  subscript  1  denotes  plate  1. 
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Figure  10.  Double  Plate  with  No  Fluid  Load 
Combining  equations  (109)  -  (116)  yields  the  eight-by-eight  linear  system  of  equations 


Ax  =  b , 

(117) 

where  the  entries  of  equation  (1 17)  are 

A,  j  =  (-a?A2  -2a\n2  -A2£^)exp(ia2c), 

(118) 

A, 2  =  (-«?A2  -2a\n2  -  A2£^)exp(-ia2c) , 

(119) 

A,3  =  (-  2kxPiF2  )exp(i/?2c) , 

(120) 

A, 4  =  . 

(121) 

A  ,5  =  0  , 

(122) 

> 

on 

II 

O 

(123) 
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4,7=0,  (124) 

4,8  =  0  »  (125) 

A2l=(-2fA2kxa2)exp(icc2c),  (126) 

4.2  =  (2M2kxa2  )exp(-i  a2c) ,  (127) 

4.3  =  (fhtf  ~  Ml  )exp(ij02c) ,  (128) 

4.4  =  (“2#  -  !hkl  )exp(-ij32c) ,  (129) 

4,5=°^  (130) 

4,6=0,  (131) 

4,7=0,  (132) 

4,8  =  0 .  (133) 

4.1  =  (i  a2  )exP(i  a2b) ,  (134) 

4.2  =  (-  i«2  )exp(-ia2&) ,  (135) 

4.3  =  (i*JexP(i/M0 ,  (136) 

4.4  =  (4)exPH4i)>  (137) 

4.5  =  (-  iai  )exp(i axb) ,  (138) 

4.6  =  (iai  )exp(-ia1fe) ,  (139) 

4.7  =  ('  ikx  )exP0A  b) ,  '  (140) 

4.8  =(-i4)exP(-iA^)’  (141) 

4,i  =(i4)exp(ia2fc),  (142) 
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A4,2  =(i^)eXPH«2&)’ 

(143) 

4*4.3  =(-i£2)exP(i/^)> 

(144) 

A4,4  =  (ift)eXP(-#). 

(145) 

4*4,5  =(-i^)exp(ia1&), 

(146) 

4*4,6  =(-i^)exp(-ia,Z>), 

(147) 

4*4,7  =  (iA  )exp(i/3,£) , 

(148) 

4*4,8  =  (~iA)exPHW> 

(149) 

4*5,1  =  (-«?4  -2a?iu2  -^^)exp(ia2fo) , 

(150) 

4*5,2  =  -2«?ju2  -A2^)exp(-ia2fc) , 

(151) 

4*5.3  =  (“  2M2/*2)exP(i&&)  > 

(152) 

4*5,4  =  (2AA^2)exPHM , 

(153) 

4*5,5  =  («?A]  +  2a?/*!  +  A1fc2)exp(ia1fc) , 

(154) 

4*5,5  =  (a;2  A,  +2a?/ij  +A]£2)exp(-ia]b), 

(155) 

4*5.7  =(2^A^i)exp(iA&). 

(156) 

4*5,8  =  (-  2*JCAtt)exPHA*0  > 

(157) 

4*6,1  =  (-2^¥2)exP(i«2fc)  . 

(158) 

4*6,2  =(2^2Aa2)exPH«2^)> 

(159) 

4*6,3  =  (hifi  -  ^2kx  )exp(i/326) , 

(160) 

4*6,4  =  )exp(-i  . 

(161) 

4*6.5  =(2^i^«i)exp(ia]&), 

(162) 

4,6  =  (-  2ftkxa\  )exp(-ia,6) , 

(163) 

4,7  =  (~ViP?  +  ftkx  )exp(i^j6) , 

(164) 

4,8  =  (“  ftP ?  +  )exp(-iA4 , 

(165) 

> 

II 

O 

(166) 

o' 

II 

<N 

(167) 

4,3  -  0  > 

(168) 

4,4  =  0, 

(169) 

4,5  =  (ra?^ i  _2a?Ju1  -Aj^  )exp(ia1(3), 

(170) 

4,6  =  (-  -  2a?Mi  -  44?  )exp(-ia,a) , 

(171) 

4,7  =  (-  2kxPift  )exp(i/3,a) , 

(172) 

4,8  =  foxPift  )exp(-iAa) » 

(173) 

II 

O 

(174) 

o 

II 

04 

(175) 

^8,3  “  0 » 

(176) 

> 

II 

O 

(177) 

4,5  =  (“  2ftkxa\  )exp(i  c^a) , 

(178) 

4,6  =  (2ftkxa\  )exp(-ia1a) , 

(179) 

4,7  =  (ftA2  -  ftkx  )exp(i Pxa) , 

(ISO) 

4,8  =  (ftP?~  ftkx  )exp(-iAa) . 

(181) 
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*u  =A2{kx,(o), 

(182) 

*2,1  =  B2(kx,(o), 

(183) 

*3.1  =  ^2  (%x  > ®)  » 

(184) 

*4,1  =  ^2  (^jc  *  » 

(185) 

*5,1  =Mkx,(0), 

(186) 

*6,1  = 

(187) 

*7,1  =  C, (kx,a>) , 

(188) 

*8,1  =  Dl(kx,0)) , 

(189) 

^1,1  —  — (^Jt  ’  ®)  > 

(190) 

^2,1  =  0  * 

(191) 

o 

II 

(192) 

o 

II 

(193) 

la 

II 

O 

(194) 

o 

II 

(195) 

o 

II 

(196) 

and 

o 

II 

oo 

(197) 

Using  equations  (1 18)  -  (197),  the  solution  to  the  constants  A, , 
can  be  found  by 

,  C„  Dj ,  Aj ,  B2 ,  C2 ,  and  D 

x  =  A_1b. 

(198) 

Additionally,  the  transfer  function  between  the  normal  displacement  at  location  z  (when  z<b) 

and  the  structural  load  can  be  written  as 
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T(kx,(o) 


Uz(z,kx,(o) 

Ps(kx,(o) 


=  Al(kx,(o)ial  exp(i  cqz) - Bx (kx,a>)ial  expHcqzH 


C1(kx, co)ikx  exp(ijSj z)  +  Dl(kx, co)ikx  exp(-i/3, z),  (199) 

or  (when  z  >  b) 


T{kx,(0 ) 


Ps(kx,co ) 


A2(kx,(o)ia2  exp(ia2z)-  B2  (kx,(o)ia2  exp(-ia2z)  + 


C2 (kx,co)ikx  exp(i fi2z)  +  D2 (kx,(o)ikx exp(-i /32z).  (200) 

The  transfer  function  between  the  normal  stress  at  location  z  (when  z<b)  and  the  structural  load 
can  be  written  as 


T (kx , m)  =  ~^:~x ^  =  A  (kx , o))(- a?Xx  -2a? Hi -A1fc2)exp(ia1z)  + 

Bl  (kx , ct>)(-  a? Aj  -  2a? Hi  ~\k?  )exp(-i ai z)  + 

C,  (kx ,  (o)(-2kx  fa  Hi )  exp(i^j  z)  +  Dx  (kx ,  (o)(2kx  fixH\)  exp(-i)81  z), 

or  (when  z  >  b) 

T(K .®)  =  T°p^x'^  -  A  -  lain, -X%k\  )exp(ia2j)  + 

B2 (kx ,G))(- a? A2  -  2a? Hi  ~ k2k? )exp(-ia2z)  + 

C2  (kx,a))(-2kxp2H2 )  exp(i  fi2z)  +  D2  (kx ,  c o)(2kx  f32H2 )  exp(-i  fi2z). 

Figure  1 1  is  a  plot  of  displacement  in  the  ^-direction  (normal)  divided  by  normal  pressure 
versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  -  c.  In  this 
example,  the  sum  of  the  thickness  of  the  two  plates  is  very  small  (h  =  0.01  m).  The  solid  line 
represents  the  thick  plate  theory  using  two  plates  (equation  (200))  and  the  x’s  represent  the  thin 
plate  theory  (equation  (68))  using  h  =  0.01  m.  Additionally,  the  material  properties  of  the  two 
(thick)  plates  are  identical  so  that  a  direct  comparison  between  the  thick  and  thin  plate  theories 
can  be  made  using  only  a  single  thin  plate.  Parameters  used  to  formulate  this  model  are:  bottom 
of  plate  1  (a)  is  -0.010  m,  intersection  of  plates  2  and  1  ( b )  is  -0.005  m,  top  of  plate  2  (c)  is  0  m. 

Young’s  modulus  ( E )  is  13.2e9  (l+0.03i)  N/m2,  Poisson’s  ratio  (v)  is  0.30,  and  density  of  the 
plate  (p)  is  1938  kg/m3 .  Only  positive  values  of  wavenumber  are  shown,  as  the  function  is 
symmetric  about  k  =  0.  Note  that  there  is  almost  complete  agreement  between  the  thick  plate 
theory  using  two  similar  thick  plates  and  the  thin  plate  theory  using  one  plate  for  this  specific 
example.  Figure  12  is  a  plot  of  displacement  in  the  z-direction  (normal)  divided  by  normal 
pressure  versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  = 
c.  In  this  example,  the  sum  of  the  thickness  of  the  two  plates  is  very  small  ( h  =  0.015  m).  The 
solid  line  represents  the  thick  plate  theory  using  two  plates  (equation  (200))  and  the  x’s  represent 


(201) 


(202) 
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the  thin  plate  theory  (equation  (68))  using  h  =  0.005  m.  Additionally,  the  material  properties  of 
the  two  (thick)  plates  are  extremely  dissimilar  so  that  a  direct  comparison  between  the  thick  and 
thin  plate  theories  can  be  made  using  only  a  single  thin  plate  with  the  thickness  of  the  stiff  thick 
plate  ( h  =  0.005  m).  Parameters  used  to  formulate  this  model  are:  bottom  of  plate  1  (a)  is  -0.015 
m,  intersection  of  plates  2  and  1  (b)  is  -0.005  m,  top  of  plate  2  (c)  is  0  m,  Young’s  modulus  of 

plate  l(£j)is  13.2e9  (l+0.03i)  N/m2,  Young’s  modulus  of  plate  2  (E2)  is  13.2e4  N/m2, 

Poisson’s  ratio  of  plates  1  and  2  (u)  is  0.30,  and  density  of  plates  1  and  2  (p)  is  1938  kg/m3 . 
Only  positive  values  of  wavenumber  are  shown,  as  the  function  is  symmetric  about  k  =  0.  Note 
that  there  is  almost  complete  agreement  between  the  thick  plate  theory  using  two  dissimilar  thick 
plates  and  the  thin  plate  theory  using  one  plate  (matched  to  the  stiff  plate  parameters)  for  this 
specific  example.  Figure  13  is  a  plot  of  the  wavenumber-frequency  (koo)  surface  of  the  system 
using  the  thick  plate  equations  of  motion  (equation  (200))  evaluated  at  z  =  c.  The  first  plate  is  a 
stiff  material  and  the  parameters  of  this  plate  are:  bottom  of  the  plate  (a)  is  -0.06  m,  top  of  the 

plate  (b)  is  -0.05  m,  Young’s  modulus  of  the  plate  (£j)  is  4.55el0  (l+0.03i)  N/m2,  Poisson’s 

ratio  of  the  plate  (u, )  is  0.30,  and  density  of  the  plate  ( px )  is  7700  kg/m3 .  The  second  plate  is  a 
soft  material  and  the  parameters  of  this  plate  are:  bottom  of  the  plate  (b)  is  -0.05  m,  top  of  the 
plate  (c)  is  0  m,  Young’s  modulus  of  the  plate  (£2)  is  le9(l+0.15i)  N/m2,  Poisson’s  ratio  of  the 

plate  (v2 )  is  0.45,  and  density  of  the  plate  ( p2 )  is  1200  kg/m3 .  The  color  scale  to  the  right  of 

the  plot  is  in  decibels.  All  values  greater  than  -190  dB  are  displayed  as  -190  dB  and  all  values 
less  than  -220  dB  are  displayed  as  -220  dB. 
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Phase  Angle  (degrees)  Magnitude  (dB) 


Figure  11.  Transfer  Function  of  a  Double  Plate  with  No  Fluid  Load 
( Equation  (200))  Using  Similar  Material  Properties  at  500  Hz 


Phase  Angle  (degrees)  Magnitude  (dB) 


Figure  12.  Transfer  Function  of  a  Double  Plate  with  No  Fluid  Load 
(Equation  (200))  Using  Dissimilar  Material  Properties  at  500  Hz 
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Frequency  (Hz) 


Figure  13.  Transfer  Surface  of  a  Double  Plate  with  No  Fluid  Load 
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6.  DOUBLE  PLATE  WITH  A  FLUID  LOAD 


The  fourth  transfer  function  derived  is  a  double  plate  with  a  fluid  load.  The  applied  load 
is  modeled  as  an  incident  pressure  wave  in  the  fluid  at  definite  wavenumber  and  frequency  at 
location  z  =  c,  as  shown  in  figure  14.  The  normal  and  tangential  stresses  in  the  system  at  the 
boundary  z  =  c  are 


r=(,,c,0  =  (A2+2ft)^%^+A2*!4M  =  -p„(x>c,0. 


and 


rzx(x,c,t)  =  fi2 


dz 


dux(x,c,t )  du  A  x,c,t) 


dz 


-+- 


dx 


dx 


=  0. 


(203) 


(204) 


where  the  subscript  2  denotes  plate  2.  The  interface  between  plates  2  and  1  requires  four 
equations.  The  first  two  are  displacement  constraints,  which  are 


plate  1  5 


and 


Ux(XMLte2=Ux(XMl 


Iplate  2  7  "'I  plate  1 

The  second  two  are  stress  constraints,  which  are 
Tzz(xMLtt2=rzz(x,b,t)\ 


(205) 

(206) 


l  plate  1  ’ 


and 


(2°7) 

M*’Mp,ale2  =  atel  *  (208) 

Finally,  the  normal  and  tangential  stresses  in  the  system  at  the  boundary  z  -  a  are 

(209) 


ratea.OMA1  +  2ft)^4^+A1^%2^  =  0, 


and 


rcc(x,a,0  =  Mi 


dz 

dux(x,a,t)  du  Ax,  a,  t) 


dz 


dx 


dx 


=  0. 


(210) 


where  the  subscript  1  denotes  plate  1. 
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Fluid 


Figure  14.  Double  Plate  with  Fluid  Load 

Combining  equations  (75)  -  (78)  and  (203)  -  (210)  yield  the  eight-by-eight  linear  system 
of  equations 

Ax  =  b,  (211) 


where  the  entries  of  equation  (211)  are 


A,i  = 


2„. 


■a\X2  -2 a\\x2  ~X2kl  +  PfC°  °2 

y 


exp(i  a2c), 


(212) 


A,2  ~ 


-a\X2  -la2p2  -X2k2x 


exp(-ia2c). 


(213) 


A,3  - 


Xkx(32fi2  + 


Pf(o2kx  ^ 


exp(i  p2c) , 


(214) 
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Pf(0  kr 

A, 4"  2M2A*2+  f  exp(  1  P2c), 

{  Y  J 

(215) 

> 

II 

o 

(216) 

A, 6  -  ^ » 

(217) 

> 

II 

O 

(218) 

A, 8 _  o , 

(219) 

A,i  =  (-  ^2kxa2  )exp(i a2c) , 

(220) 

A, 2  =  (2fA2kxa2)ex-p(-ia2c) , 

(221) 

A, 3  =  (fhP?  -  V 2kx  )exp(i^2c)  , 

(222) 

A, 4  =  ifhfi  -  ^2kl  )exp(-iy32c) , 

(223) 

A,5  =  0  > 

(224) 

A, 6 =  o  ’ 

(225) 

A,7  =  0 , 

(226) 

o 

II 

OO 

(227) 

A,i  =  (i«2)exp(ia2&) , 

(228) 

A,2  =  (-ia2)exp(— ia2b). 

(229) 

A, 3  =  (ikx)exv(ip2b) , 

(230) 

A,4  =  (iA)exPHA&), 

(231) 

A  ,5  =(-ia1)exp(ia1ft). 

(232) 

A, 6  =(ia1)exp(-ia1fe), 

(233) 
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4,7  =(-i^)exp(i^), 

(234) 

4,8  =(-ikx)exp(-iPlb), 

(235) 

4,i  =  {ikx)exip(ia2b) , 

(236) 

4,2  =  (i^)exp(-ia2fe) , 

(237) 

4,3  =(-i4)exp(i4^), 

(238) 

4,4  =  (i/?2)exp(-ij326) . 

(239) 

4,5  =(-i^)exp(ia1fe), 

(240) 

4,6  =(-ikx)exp(-ia1b), 

(241) 

4,7  =  (i  A  )exp(i/3,&) , 

(242) 

4,8  =  (-iA)expHA^)> 

(243) 

4,i  =  (-a?4  -2a?/i2  -  4**  )exp(ia2&) , 

(244) 

4,2  =  (-afAj  -2alpL2  ~4^)exp(-ia2t) , 

(245) 

4,3  =  (- 24^2^2  )exP(i A>&), 

(246) 

4,4  =  (2^02^2  )exP(-i&6) , 

(247) 

4,5  =  +  2a?^  +  Aj^jexpGa^) , 

(248) 

4,6  =  +  2a?/z,  +  Aj^jexpHa,^) , 

(249) 

4,7  =(2^A^i)exp(iAfe). 

(250) 

4,8  =  (~  24  A  Mi  )exp(-i/3,^) , 

(251) 

4,i  =  ("  ^2kxa2  )exp(ia2fe) , 

(252) 
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4.2  =  (2P 2kxa2  )exP(-i  a2b)  >  (253) 

4.3  =  -  P-2 kl )exp(i£,6) ,  (254) 

4.4  =  iPlPz  ~  p2kx  )®XP(— i&2b)  ’  (255) 

4.5  =  (2Pikxa\  )exP(i  a\b)  >  (256) 

4.6  =  (“  2P\kxa\  )exPH  «i  b) ,  (257) 

4.7  =  (“  P\P'  +  Pikx  )exp(i p{b) ,  (258) 

4.8  =  (“  PiP?  +  Pxkl  )exp(-iA&) ,  (259) 

4,i=  0^  (260) 

4,2=0,  (261) 

4.3  =  0,  (262) 

4.4  =  0,  .  (263) 

4.5  =  (-  «?4  -  2a?/U]  -  )exp(iaja) ,  (264) 

4.6  =  (-  -  2ai>,  -  4^  )exp(-ia,a) ,  (265) 

4.7  =  (~  2kxP\Pi  )exp(i/31fl) ,  (266) 

4.8  =  (2kxP\Pi  )exp(-iy31a) ,  (267) 

4,i=  °>  (268) 

4,?=°’  ,  (269) 

4,3  =  0,  (270) 

4,4=0,  (271) 

4,5  =  (-2^!^!  )exp(i  OTja),  (272) 
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As.6  =  {2ftkxai  JexpHOjfl) , 


(273) 


=  (M?  -^A2)exp(i^fl), 

(274) 

As,8  =  (^l^1  ~  )exP(— lP\a)  > 

(275) 

*1,1  =A2(kx,0)), 

(276) 

*2,1  =B2(kx,oj), 

(277) 

*3,1  =  ^2  (^jt  > j 

(278) 

*4,1  =D2(kx,(0), 

(279) 

*5,1  =Al(kx,(0), 

(280) 

*6,1  ='B1(^,W), 

(281) 

*7,1  =Ci(*jt,  ®). 

(282) 

*8,1  =  A  (^,W), 

(283) 

*>u  =  -27/  (£*,«)  exp(-iyc) , 

(284) 

K) 

II 

o 

(285) 

o 

II 

(286) 

jf4 

II 

o 

(287) 

cy- 

<Jr 

II 

o 

(288) 

o 

II 

vo 

-Cj 

(289) 

O 

II 

(290) 

oo 

II 

O 

(291) 
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Using  equations  (212)  -  (291),  the  solution  to  the  constants  Ax,  Bx,  Cp  Dx,  A2,  B2,  C2,  and  D2 
can  be  found  by 

x  =  A-1b .  (292) 

Additionally,  the  transfer  function  between  the  normal  displacement  at  location  z  (when  z  <  b) 
and  the  structural  load  can  be  written  as 


nK '“) = Zm = A  °ma'z)  " 

Bx  (kx ,  co)i  ax  exp(-i  ax  z)  + 

C1(kx,co)ikx  exp(ij3,z)  +  Dl(kx,co)ikx  exp(-i/3,z),  (293) 

or  (when  z>b) 

T ( kx,o) )  =  — ^ ,<y)ia2  exp(ia2z)  - 
Pi(kx,o>)  exp(-iyc) 

B2{kx,(o)\a2  exp(-i  a2z)  + 

< 

C2(kx,o))ikx  exp(i/32z)  +  D2 (kx,o))ikx  exp(-i/l2z).  (294) 


The  transfer  function  between  the  normal  stress  at  location  z  (when  z<b)  and  the  structural  load 
can  be  written  as 

T  (kx  ,co)  =  =  A2  (kx,co)ia2  exp(i  a2z)  - 

P,(kx,a))exv(-\Yc) 

B2(kx,co)ia2cxp(-ia2z)  + 

C2 (kx,(o)ikx  exp(i fi2z)  +  D2 (kx,co)ikx  exp(-i/J2z).  (295) 


or  (when  z  >  b) 


T(kx,co)  =  - (^^^p^iyc)  =  A2(kx,a)){-a22X2  -2 a\ix2  -A2fc;)exp(ia2z)  + 
B2(kx,eo)(-a;X2  -2a2n2  -k2k2x  )exp(-i a2z)  + 

C2  (kx ,  (o){-2kxP2jj,2 )  exp(i/32  z)  +  D2  (kx  ,eo)(2kxp2n2 )  exp(-i£2z).  (296) 
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Figure  15  is  a  plot  of  the  displacement  in  the  z-direction  (normal)  divided  by  normal 
pressure  versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  = 
c.  In  this  example,  the  sum  of  the  thickness  of  the  two  plates  is  very  small  ( h  =  0.01  m).  The 
solid  line  represents  the  thick  plate  theory  using  two  plates  (equation  (294))  and  the  x’s  represent 
the  thin  plate  theory  (equation  (107))  using  h  =  0.01  m.  Additionally,  the  material  properties  of 
the  two  (thick)  plates  are  identical  so  that  a  direct  comparison  between  the  thick  and  thin  plate 
theories  can  be  made  using  only  a  single  thin  plate.  Parameters  used  to  formulate  this  model  are: 
bottom  of  plate  1  (a)  is  -0.010  m,  intersection  of  plates  2  and  1  (b)  is  -0.005  m,  top  of  plate  2  (c) 

is  0  m,  Young’s  modulus  (£)  is  13.2e9  (l+0.03i)  N/m2,  Poisson’s  ratio  (v)  is  0.30,  and  density 
of  the  plate  (p)  is  1938  kg/m3 .  Only  positive  values  of  wavenumber  are  shown,  as  the  function 
is  symmetric  about  k  =  0.  Note  that  there  is  almost  complete  agreement  between  the  thick  plate 
theory  using  two  similar  thick  plates  and  the  thin  plate  theory  using  one  plate  for  this  specific 
example.  Figure  16  is  a  plot  of  the  displacement  in  the  z-direction  (normal)  divided  by  normal 
pressure  versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  = 
c.  In  this  example,  the  sum  of  the  thickness  of  the  two  plates  is  very  small  ( h  =  0.015  m).  The 
solid  line  represents  the  thick  plate  theory  using  two  plates  (equation  (294))  and  the  x’s  represent 
the  thin  plate  theory  (equation  (107))  using  h  =  0.005  m.  Additionally,  the  material  properties  of 
the  two  (thick)  plates  are  extremely  dissimilar  so  that  a  direct  comparison  between  the  thick  and 
thin  plate  theories  can  be  made  using  only  a  single  thin  plate  with  the  thickness  of  the  stiff  thick 
plate  ( h  =  0.005  m).  Parameters  used  to  formulate  this  model  are:  bottom  of  plate  1  (a)  is  -0.015 
m,  intersection  of  plates  2  and  1  (b)  is  -0.005  m,  top  of  plate  2  (c)  is  0  m,  Young’s  modulus  of 

plate  2  (E2)  is  13.2e9  (l+0.03i)  N/m2,  Young’s  modulus  of  plate  1  (Ej)  is  13.2e4  N/m2, 

Poisson’s  ratio  of  plates  1  and  2  ( v )  is  0.30,  and  density  of  plates  1  and  2  (p)  is  1938  kg/m3 . 

Only  positive  values  of  wavenumber  are  shown,  as  the  function  is  symmetric  about  k  =  0.  Note 
that  there  is  almost  complete  agreement  between  the  thick  plate  theory  using  two  dissimilar  thick 
plates  and  the  thin  plate  theory  using  one  plate  (matched  to  the  stiff  plate  parameters)  for  this 
specific  example.  Figure  17  is  a  plot  of  the  wavenumber-frequency  (kco)  surface  of  the  system 
using  the  thick  plate  equations  of  motion  (equation  (294))  evaluated  at  z  =  c.  The  first  plate  is  a 
stiff  material  and  the  parameters  of  this  plate  are:  bottom  of  the  plate  (a)  is  -0.06  m,  top  of  the 

plate  ( b )  is  -0.05  m,  Young’s  modulus  of  the  plate  (E,)  is  4.55el0  (l+0.03i)  N/m2,  Poisson’s 

ratio  of  the  plate  (vt )  is  0.30,  and  density  of  the  plate  ( p, )  is  7700  kg/m3 .  The  second  plate  is  a 
soft  material  and  the  parameters  of  this  plate  are:  bottom  of  the  plate  (b)  is  -0.05  m,  top  of  the 
plate  (c)  is  0  m,  Young’s  modulus  of  the  plate  (E2)  is  le9  (l+0.15i)  N/m2,  Poisson’s  ratio  of  the 

plate  (v2 )  is  0.45,  and  density  of  the  plate  ( p2 )  is  1200  kg/m3 .  The  color  scale  to  the  right  of 

the  plot  is  in  decibels.  All  values  greater  than  -190  dB  are  displayed  as  -190  dB  and  all  values 
less  than  -220  dB  are  displayed  as  -220  dB. 


42 


Phase  Angle  (degrees)  Magnitude  (dB) 


Figure  15.  Transfer  Function  of  a  Double  Plate  with  Fluid  Load 
(Equation  (294))  Using  Similar  Material  Properties  at  500  Hz 
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Figure  16.  Transfer  Function  of  a  Double  Plate  with  Fluid  Load 
(Equation  (294))  Using  Dissimilar  Material  Properties  at  500  Hz 
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Figure  1 7.  Transfer  Surface  of  a  Double  Plate  with  Fluid  Load 
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7.  TRIPLE  PLATE  WITH  NO  FLUID  LOAD 


The  fifth  transfer  function  derived  is  a  triple  plate  with  no  fluid  load.  The  structural  load 
consists  of  an  applied  pressure  at  definite  wavenumber  and  frequency  at  location  z  =  d,  as  shown 
in  figure  18.  The  normal  and  tangential  stresses  in  the  system  at  the  boundary  z  =  d  are 


*=  (*.  d,ty =  +  2ih)^f’t)+  A,  =  -Ps(x,d,t), 

<JZ  OX 


and 


rzx(x,d,t)  =  fi 3 


duj^dd)  ^  du^jx^dd) 
dz  dx 


=  0. 


(297) 


(298) 


where  the  subscript  3  denotes  plate  3.  The  interface  between  plates  3  and  2  requires  four 
equations.  The  first  two  are  displacement  constraints,  which  are 

UZ^CAVX ate3  =  Opiate  2  ’  (299) 

and 

=  ■  (300) 


The  second  two  are  stress  constraints,  which  are 

Tzz(x’c’ 0|plate3  =  7K(^,c,f)|plate2  ,  (301) 

and 

T«(^^0|plate3=Tzx(x,c,0|plate2  .  (302) 

The  interface  between  plates  2  and  1  requires  four  equations.  The  first  two  are  displacement 
constraints,  which  are 

u^x’b’,\^l=u^x’b’,\^  .  (303) 

and 

“^’Mplate2  =  “*(*’MpIatei  •  (304) 


The  second  two  are  stress  constraints,  which  are 


plate  1 


(305) 


and 
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Figure  18.  Triple  Plate  with  No  Fluid  Load 


M*>Mplate2=M*>Mplate ,  • 

Finally,  the  normal  and  tangential  stresses  in  the  system  at  the  boundary  z  =  a  are 


(306) 


and 


tzx(x,a,t)  =  flx 


dz  dx 


dux(x,a,t)  |  <faz(x,a,Q 


(307) 

(308) 


where  the  subscript  1  denotes  plate  1. 
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Combining  equations  (297)  -  (308)  yields  the  twelve-by-twelve  linear  system  of  equations 
Ax  =  b,  (309) 

where  the  entries  of  equation  (309)  are 

Au  =  (-  a32  A 3  -  2a2fi3  -  A3 kl  )exp(i  a3d) ,  (310) 

A12  =  (-o^Aj  -  2a} fi3  ~^k2 )exp(-ia3rf),  (311) 

A,3=(-2^^3^3)exP  (312) 

AlA  =  {2kJ3n3)exV(-ip3d),  (313) 

4,5=0,  (314) 

4,6=0,  (315) 

4,7=0,  (316) 

4,8  =°,  •  (317) 

A19=0,  (318) 

Ayo=0,  (319) 

Ain  =  0,  (320) 

AU2  =  0 ,  (321) 

4.1  =  (“  )exp(i  a3d) ,  (322) 

4.2  =  (2^34^3  )exp(-ia3rf) ,  (323) 

A2^(jx3fi}-/i3kl)tx^3d),  (324) 

4,4  =  ^/3i-fx3k2x)exV(-il33d),  (325) 

A2>5  =0,  (326) 
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4.6  ~  0 , 

(327) 

4,7  =  0 , 

(328) 

cT 

II 

oo 

(329) 

JO 

VO 

II 

o 

(330) 

4,10  “0» 

(331) 

JO 

II 

O 

(332) 

4,12  ~  0 , 

(333) 

4,i  =  (ia3)exp(ia3c) , 

(334) 

4,2  =(-ia3)exp(-ia3c). 

(335) 

4s  =  (ifcJexpCi&c) , 

(336) 

4,4  =  (i4)exp(-i/33c), 

(337) 

4,5  =  (-ia2)exp(ia2c) , 

(338) 

4,6  =(i«2)exp(-i  a2c). 

(339) 

4,7  =(-i4)exP0^2c). 

(340) 

4,8=(-i4)exPH^2c), 

(341) 

4,9  =  0  > 

(342) 

4,10 =  4 

(343) 

4,ii =  4 

(344) 

4,12  =  o , 

(345) 

4,1  =  (i£x)exp(ia3c), 

(346) 
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4,2  =  k)eXP("i«3C)> 

(347) 

A4,3=(-i^3)eXP(i^3C). 

(348) 

Am  =  (i£3)exp(-i/?3c), 

(349) 

A4,5  =(-i^)eXP(i«2C)- 

(350) 

A4,6  =  (-i**  )exp(-i  a2c), 

(351) 

A4,7  =  (i/^2  )exP(i/^2C)  » 

(352) 

A4,8  ~  (_  i/^2  )exp(-i/32c)  , 

(353) 

A*,9  =  0  > 

(354) 

o 

1! 

o 

(355) 

£ 

II 

o 

(356) 

o 

II 

<N 

(357) 

A;.,  =  (-a32A 3  -  2a32^3  -/I3/:2  )exp(i  a3c), 

(358) 

^5,2  =  (-«3%  -  2a 2  ^3  —  ^3^2  )exp(-i  a3c) , 

(359) 

A.,3  =  (-  2kxP:^3  )exp(i/V) . 

(360) 

A,4  =  (2^/33iu3)exp(-i^3c) , 

(361) 

^5,5  =  («?4  +  2a?/42  +  )exp(i  a2c) , 

(362) 

A, 6  =  (fxi*2  +  2a^2  +  ^2kx  )exp(-i  a2c) , 

(363) 

A>,7  =  (2£*/?2jU2)exp(ij32c), 

(364) 

^5,8  =  (“  2kxPi^2  )exp(-ij32c) , 

(365) 

2^5  9  =  0  , 

(366) 
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4,io  “O’ 


(367) 


4,n  “  0 » 

(368) 

K> 

II 

o 

(369) 

Ay  =  (-  2^3kxa3  )exp(ia3c) , 

(370) 

A, 2  =  {2/^3kxa3)e\p(-ia3c) , 

(371) 

4,3  =  fa#  -  A*342  )exp(i/33c) , 

(372) 

4,4  =  (M2  "  )exp(-iy33c) , 

(373) 

4,5  =  (27*2^2  )exP(ia2c)  ’ 

(374) 

4,6  =  (“  2  AJ2  kxa2 )  exp(-i  a2c) , 

(375) 

4,7  =  (“  ^20 2  +  7*2  kx  )exp(i/32c) , 

(376) 

4,8  =  (“  2  +  f^2kx  )exp(-ijS2c) , 

(377) 

^6,9  =  ®  » 

(378) 

O 

II 

o 

(379) 

> 

II 

O 

(380) 

4,12 -  o  ’ 

(381) 

II 

O 

(382) 

> 

io 

ll 

o 

(383) 

4,3  “  o , 

(384) 

> 

II 

O 

(385) 

4,5  =  (ia2)exp(io;2&), 

(386) 
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4,6  =  (-  ia2  )exp(-i  a2b) ,  (387) 

4 j  =  (ikx)exp(ifi2b) ,  (388) 

48  =  k)exP(-#).  (389) 

4.9  =  {-iax  )exp(i  axb) ,  (390) 

4.10  =  (iai  )exp(-ia,6) ,  (391) 

4.1  i  =  ('  ikx  )exp(iA&) ,  (392) 

A7>i2  =(-i^)exp(-i^),  (393) 

4.1  =0,  (394) 

4.2  =0,  (395) 

4.3  =  0,  (396) 

4.4  =  0,  (397) 

4.5  =  (4  )exp(ia2&) ,  (398) 

4.6  =  (ikx  )exp(-i  a2b) ,  (399) 

4.7  =  (-iA>)exP(iA>&)>  (400) 

4.8  =(i^2)exp(-ijS2&),  (401) 

4.9  =  ('ikx  )exp(iccxb) ,  (402) 

4.10  =  (-  ikx  )exp(-i  axb) ,  (403) 

4,i  i  =  (i  A  )exp(i/3jib) ,  (404) 

4,12  =  (“  *A  )exp(-iA^) .  (405) 

4,i=0,  (406) 
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A), 2  —  0  , 

(407) 

> 

OJ 

II 

o 

(408) 

> 

V 

II 

o 

(409) 

A, 5  =  (~  ccl^  -2 a?yU2  -  A2^)exp(ia2^)  * 

(410) 

A, 6  =  -2a?^2  -A2^)exp(-ia2d) , 

(411) 

A>.7  =  (-  2kxP2n2)exp(i/32b) , 

(412) 

A,8=(2^ftft)exPHftJ)' 

(413) 

A)  ,9  =  (a?A  +  2a?/i,  +A1^2)exp(ia]&) , 

(414) 

A, 10  =  (a?A  +2a?iu1  +  A&x)exp(-iaj&), 

(415) 

Aui  =(2^Aiui)exp(i/31Z?), 

(416) 

Aw  =  (-2^AMi)exp(-i^fe), 

(417) 

> 

p 

II 

o 

(418) 

> 

p 

N> 

II 

o 

(419) 

Ao,3  =  0  * 

(420) 

Ao,4  =  ^  > 

(421) 

A  0,5  =  (-  2^A«2  )exp(ia2&) , 

(422) 

A  0,6  =(2Ai2A«2)exPH«2^)> 

(423) 

A 0,7  =  (M*  -^2^2)eXP(iA&)  . 

(424) 

A  0,8  =  if1 2  ft  ~  l*2kx)&XVHp2b')  . 

(425) 

A  0,9  =  (2^iAai)exP(iai&). 

(426) 
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Ao.io  =  {-2nlkxa1)exp(-ialb), 
Ao.il  =  (“  Mi  A2  +  Mi kl  )exp(i P}b) , 
Ao,i2  =  (-  Mi#2  +  Mi  *1  )exp(-i  #Z>) , 


A  u  =  0 , 

A  1,2 =  o  > 

A  i,3 =  0’ 

A  1,4 

A  1,5  =^’ 

A  1,6  =  0  ’ 


A  1,7  ==  0  > 

A  1,8  = 

A i,9  =  (- a.2  A  -2a,2Ju1  -\k2x  )exp(iaja), 
A  uo  =  (-a?  A  “2  a?/4j  -A^)exp(-iaja), 
Aui  =  (~  2kx  AMi  )exp(i#a) , 

A  i,i2  =  faxfiifh  )exp(-iAa)  > 

A  2,1  =  0  > 

A 2,2  =  0  > 


A  2,3  -  0  > 
A  2,4  =  0  > 


(427) 

(428) 

(429) 

(430) 

(431) 

(432) 

(433) 

(434) 

(435) 

(436) 

(437) 

(438) 

(439) 

(440) 

(441) 

(442) 

(443) 

(444) 

(445) 
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(446) 


•At  2,5  “  0  * 

A2,6=0>  (447) 

A  2,7  =  0  >  (448) 

A2,8=0’  (449) 

A  2,9  =  ("  2^\kxa\  )exp(ia1a) ,  (450) 

A240  =  (2^i^a,)exp(-ia1a),  (451) 

A 2,1 1  =  U#  ~  Vikx  )exp(i .j8,a) ,  (452) 

A2.12  =  (ftP?  -  A^^JexpHfta),  (453) 

*1,1  =  (454) 

*2,1  =B3(kx,(o),  (455) 

*3,1  =C3(kx,(0),  (456) 

*4,1  =D3(kx,(0),  (457) 

*5,1  =A2(kx,(D),  (458) 

*6,1  =B2(kx,(0),  (459) 

*7,1  =  C2  (kx » »  '  (460) 

*8,1  =  D2{kx,(0),  (461) 

*9,i  =  A  (**,<»),  (462) 

*10,1  =  A  (*,.<»),  (463) 

*11,1  =  C}(kx,a)),  (464) 

*12,1  =  (kx,co) ,  (465) 
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K\  = -ps(kx’°>)  ’ 

(466) 

o 

II 

(467) 

© 

II 

(468) 

o 

II 

(469) 

y» 

II 

O 

(470) 

&6,1=V 

(471) 

O 

II 

*cT 

(472) 

o 

II 

OO 

(473) 

II 

O 

(474) 

JS- 

© 

II 

o 

(475) 

bni  -0, 

(476) 

and 

JP1 

N> 

II 

o 

(477) 

Using  equations  (310)  -  (477),  the  solution  to  the  constants  Al ,  Bx,  Cx,  Dx,  A2,  B2, 
A3,  B3 ,  C3 ,  and  D3  can  be  found  by 

C  2  »  ^2  » 

II 

> 

1 

H- * 

c r 

(478) 

Additionally,  the  transfer  function  between  the  normal  displacement  at  location  z  (when  z<b) 
and  the  structural  load  can  be  written  as 

T(kx,(o)  =  *  *  =A1(^,a>)ia]exp(ia1z)-51(^,ft))ict1exp(-ia1z)  + 

•*5  >  &)) 

C1(kx,(o)ikx  exp(ijSjZ)  +  Dx  (kx,(o)ikx  exp(-i  pxz). 

(479) 

or  (when  z>b  and z<c) 

57 

T (kx ,(o )  =  — - ’  —  =  A2 (kx , a>)ia2  exp(i a2 z ) -  B2 (kx,co) i a2  exp(-i  a2z)  + 


Ps(kx,(o) 

C2 (kx,(o)ikx  exp(i p2z)  +  D2 (kx,(o)\kx  exp(-i/32z). 


(480) 


or  (when  z>  c) 


T(kx,(0)  =  =  A3(kx,(o)ia3  exp(i  a3z)-  B3{kx,(o)\a3  exp(-i  a3z)  + 

rs  \KX,(D) 

C3 (kx,(0)ikx  exp(i£3z)  +  D3 (kx,co)ikx  exp(-i p3z). 


(481) 


The  transfer  function  between  the  normal  stress  at  location  z  (when  z<b)  and  the  structural  load 
can  be  written  as 


T{kx,(o)  =  =  Al(kx,co)(-afA1  -2 \ccfni - \k2x)txp(ialZ)  + 

B\  (kx ,(o)(- a12Al  -  2a2/xl  - \k]  )exp(-iajz)  + 

Cx(kx,  (o)(-2kx/31  ^ )  exp(iy3j  z)  + Dl(kx,(o)(2kxp]/x1)  exp(-i£,  z). 


(482) 


or  (when  z  >  b  and  z<c) 

T ( kx,co )  =  ~ ’  ^  =  A, (kx,co)(- a\X2  -  2a2ji2  -A2£2)exp(ia,z)  + 

rs  \kx,CO) 

B2(kx,co)(-a2A.2  -2 a;n2  -X2k2x)exp{-\a2z)  + 

A  (kx ,  (0)(-2kx  /32/t2 )  exp(i/3,z)  +  D2  (kx  ,(D)(2kxp2/u2)  exp(-i  j32  z),  (483) 

or  (when  z>  c) 

nK .®)  =  ^  (*,.<»)(-«,%-  2a,v,  -  A,*;  )exp(i  a,  z)  + 

53(^x,ft))(-a3A3  -2 a3Ju3  -A3^)exp(-ia3z)  + 

C3  (kx,(o)(-2kxP3n3 )  exp(i/?3  z)  +  D3  (*x ,  (o)(2kxp3/x3 )  exp(-ijS3  z).  (484) 
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Figure  19  is  a  plot  of  the  displacement  in  the  z-direction  (normal)  divided  by  normal 
pressure  versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  = 
d.  In  this  example,  the  sum  of  the  thickness  of  the  three  plates  is  very  small  ( h  =  0.045  m).  The 
solid  line  represents  the  thick  plate  theory  using  three  plates  (equation  (484))  and  the  x’s 
represent  the  thin  plate  theory  (equation  (68))  using  h  =  0.045  m.  Additionally,  the  material 
properties  of  the  three  (thick)  plates  are  identical  so  that  a  direct  comparison  between  the  thick 
and  thin  plate  theories  can  be  made  using  only  a  single  thin  plate.  Parameters  used  to  formulate 
this  model  are:  bottom  of  plate  1  (a)  is  -0.045  m,  intersection  of  plates  2  and  1  (b)  is  -0.030  m, 
intersection  of  plates  3  and  2  ( c )  is  -0.015  m,  top  of  plate  3  (d)  is  0  m,  Young’s  modulus  (£)  is 

13.2e9  (l+0.03i)  N/m2,  Poisson’s  ratio  (t>)  is  0.30,  and  density  of  the  plate  (p)  is  1938  kg/m3 . 
Only  positive  values  of  wavenumber  are  shown,  as  the  function  is  symmetric  about  k  =  0.  Note 
that  there  is  almost  complete  agreement  between  the  thick  plate  theory  using  three  similar  thick 
plates  and  the  thin  plate  theory  using  one  plate  for  this  specific  example.  Figure  20  is  a  plot  of 
the  displacement  in  the  z-direction  (normal)  divided  by  normal  pressure  versus  wavenumber  in 
the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  =  d.  In  this  example,  the  sum  of 
the  thickness  of  the  three  plates  is  very  small  ( h  =  0.0032  m).  The  solid  line  represents  the  thick 
plate  theory  using  three  plates  (equation  (484))  and  the  x’s  represent  the  thin  plate  theory 
(equation  (68))  using  h  =  0.0030  m.  Additionally,  the  material  properties  of  the  three  (thick) 
plates  are  extremely  dissimilar  so  that  a  direct  comparison  between  the  thick  and  thin  plate 
theories  can  be  made  using  only  a  single  thin  plate  with  the  thickness  of  the  stiff  thick  plate  ( h  = 
0.0030  m).  Parameters  used  to  formulate  this  model  are:  bottom  of  plate  1(a)  is  -0.0032  m, 
intersection  of  plates  2  and  1  (b)  is  -0.0031  m,  intersection  of  plates  3  and  2  (c)  is  -0.0030  m,  top 

of  plate  3  (d)  is  0  m,  Young’s  modulus  of  plate  3  (£3)  is  13.2e9  (1+0.031)  N/m2,  Young’s 
modulus  of  plates  2  and  1  (E2  and  £,)  is  13.2e4  N/m2,  Poisson’s  ratio  of  plates  1, 2,  and  3  (v) 

is  0.30;  and  density  of  plates  1,  2,  and  3  (p)  is  1938  kg/m3 .  Only  positive  values  of  wavenumber 
are  shown,  as  the  function  is  symmetric  about  k  =  0.  Note  that  there  is  almost  complete 
agreement  between  the  thick  plate  theory  using  three  dissimilar  thick  plates  and  the  thin  plate 
theory  using  one  plate  (matched  to  the  stiff  plate  parameters)  for  this  specific  example.  Figure 
21  is  a  plot  of  the  wavenumber-frequency  (kco)  surface  of  the  system  using  the  thick  plate 
equations  of  motion  (equation  (484))  evaluated  at  z  =  d.  The  first  plate  is  a  stiff  material  and  the 
parameters  of  this  plate  are:  bottom  of  the  plate  (a)  is  -0.08  m,  top  of  the  plate  (b)  is  -0.06  m, 

Young’s  modulus  of  the  plate  (£j)  is  4.55el0  (l+0.03i)  N/m2,  Poisson’s  ratio  of  the  plate  (t>, ) 

is  0.30,  and  density  of  the  plate  (pt)  is  7700  kg/m3 .  The  second  plate  is  a  soft  material  and  the 
parameters  of  this  plate  are:  bottom  of  the  plate  (b)  is  -0.06  m,  top  of  the  plate  (c)  is  -0.01  m, 
Young’s  modulus  of  the  plate  ( E2 )  is  le9(l+0. 15i)  N/m2 ,  Poisson’s  ratio  of  the  plate  (v2 )  is 

0.45,  and  density  of  the  plate  ( p2 )  is  1200  kg/m3 .  The  third  plate  is  a  stiff  material  and  the 
parameters  of  this  plate  are:  bottom  of  the  plate  (c)  is  -0.01  m,  top  of  the  plate  (d)  is  0  m, 

Young’s  modulus  of  the  plate  (£3)  is  13.2e9  (l+0.03i)  N/m2,  Poisson’s  ratio  of  the  plate  (i>3)  is 

0.30,  and  density  of  the  plate  ( p3)  is  1938  kg/m3 .  The  color  scale  to  the  right  of  the  plot  is  in 

decibels.  All  values  greater  than  -190  dB  are  displayed  as  -190  dB  and  all  values  less  than  -220 
dB  are  displayed  as  -220  dB. 


59 


Wavenumber  (rad/m) 


Figure  20.  Transfer  Function  of  a  Triple  Plate  with  No  Fluid  Load 
(Equation  (481))  Using  Dissimilar  Material  Properties  at  500  Hz 


Frequency  (Hz) 
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8.  TRIPLE  PLATE  WITH  A  FLUID  LOAD 


The  sixth  transfer  function  derived  is  a  triple  plate  with  a  fluid  load.  The  applied  load  is 
modeled  as  an  incident  pressure  wave  in  the  fluid  at  definite  wavenumber  and  frequency  at 
location  z  =  d,  as  shown  in  figure  22.  The  normal  and  tangential  stresses  in  the  system  at  the 
boundary  z  =  d  are 


and 


rzx(x,d,t)  = 


dz 


dux(x,d,t )  cfaz(x,<j,f) 


dz 


dx 


dx 


=  0, 


(485) 


(486) 


where  the  subscript  3  denotes  plate  3.  The  interface  between  plates  3  and  2  requires  four 
equations.  The  first  two  are  displacement  constraints,  which  are 


M^’C’')|plate3=Mz(*’C4 


plate  2  9 


and 


^’C’f)|plate3=M^X’C’?)| 


plate  2 


(487) 

(488) 


The  second  two  are  stress  constraints,  which  are 

T*(W)|plate3  =  MW)|plate2  ,  (489) 

and 

Tz*(^,0|plate3  =  Mw)|plate2  •  (490) 

The  interface  between  plates  2  and  1  requires  four  equations.  The  first  two  are  displacement 
constraints,  which  are 

MZ(X’Mpla«e2=“^Mpla<el  ’  (491> 

and 

“*(X’Mplate2  =  M*(X’^)|platel  ‘  (492> 


The  second  two  are  stress  constraints,  which  are 


MX’Mp,ate2  =Tzz(X’b>0l 


plate  1  ’ 


(493) 
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Figure  22.  Triple  Plate  with  Fluid  Load 
and 

r„(j:,i.,r)|pUie2  =  -r„(^,fc,r)|p]ite]  . 

Finally,  the  normal  and  tangential  stresses  in  the  system  at  the  boundary 


and 


Ta(*,a,o =(^ +2ft)ftz(;*;a,r)  =0> 

etc 

,  .  [ <&*(>, a, r)  du  7  (x,  iz,  r) 

=  W - - - +  -^T -  =0, 

ofc  cfx  _ 


where  the  subscript  1  denotes  plate  1. 


Combining  equations  (75)  -  (78)  and  (485)  -  (496)  yields  the  twelve-by-twelve  linear  system  of 
equations 

Ax  =  b,  (497) 

where  the  entries  of  equation  (497)  are 


,,  ,  ,  ,  PfCQ a3 

Ai=  -  «3  A  -  2affi3  -  ?ukt  +  -4 - -  exp(ia3  d) , 

Y 


,  ,  ,  pf(o  a-, 

A  2  =  ~«3  A  -  2afp3  -  Aj*;  — - - 1  exp(-ia3rf) , 

7 


(  pfo)2kx  ] 

A,3  =  -  2kxPi^3  +  — -  exp(i P3d) , 

r 


(498) 

(499) 

(500) 


„  p  f  co kx 

A,4  =  2kx03fi3+— - 1  exp  (-ifcd),  (501) 

{  y  j 

A,s=0>  (502) 

A,6=05  (503) 

A,7=0>  (504) 

A,8  =  0 »  (505) 

A,9=0>  (506) 

A,io  =  0 .  (507) 

A,n=0’  (508) 

A,i2  =0.  (509) 

A,i  =  (“  2p3kxa3)exp(ia3d) ,  (510) 

A, 2  =  {2t*3kxa3  )exp(-ia3d) ,  (511) 
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A>,3  =  (fhP?  ~ lhk*)exp(ip3d) , 

(512) 

aia  =  ifhP*  ~  )exp(-i P3d) , 

(513) 

> 

1! 

O 

(514) 

A2,6  ~  0  > 

(515) 

K> 

II 

o 

(516) 

<3 

II 

oo 

2^ 

^3 

(517) 

O 

II 

ON 

(518) 

A2,10  ~  0 , 

(519) 

o 

II 

^3 

(520) 

A2,12  ~  ®  ’ 

(521) 

a3,i  =  (ia3)exp(ia3c) , 

(522) 

A3,2  =(-i«3)exp(-ia3c). 

(523) 

A3.3  =  (i^)exP(i/53c), 

(524) 

A3,4  =(i**)exp(-i/53c). 

(525) 

A3  5  =  (-  i  a2  )exp(ia2c) , 

(526) 

A3,6  =(i«2)exP(-i«2c)> 

(527) 

A3,7  =(-i^)exp(ij02c). 

(528) 

A3,s  =  (“  )exp(— i/?2c) , 


(529) 


4io=  0-  •  (531) 

4n=0,  (532) 

A,  12=0,  (533) 

\i  =  (i  kx  )exp(i  a3c) ,  (534) 

A4,2  =  (ikx  )exp(-i  a3c) ,  (535) 

A4,3  =  (-  )exp(i  p3c) ,  (536) 

A4,4  =  (i^3  )exp(-i  03c) ,  (537) 

A4,5  =  (-  ikx  )exp(ia2c) ,  (538) 

A4,6  =  (-  ikx  )exp(-i  a2c) ,  (539) 

\j  =  (i& )exp(i/32c) ,  (540) 

A4,8  =  (-  i^2  )exp(-i/32c) ,  (541) 

A49  =0,  (542) 

A4j10  =0,  (543) 

A4,„=0,  (544) 

A4,12  =  0 ,  (545) 

^5,1  =  (“  ~  2«32^3  ~^3kl  )exp(i a3c) ,  (546) 

^5,2  =  (-  ~  2«32 ^3  -  hk\  )exp(-i  a3c) ,  (547) 

^5,3  =  (-  2kx^3  )exp(ij 83c)  ,  (548) 

^5,4  =  foxPiVi  )exp(-i/33c) ,  (549) 
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(550) 


4.5  =  (a  fAj  +  2  a?n2  +  A2fc^  )exp(ia2c)  > 

4.6  =  (a?4  +  2«?M2  +  ^2kl  )exp(-ia2c) ,  (551) 

4.7  =  foxPiVi  )exp(ij02c) ,  (552) 

4.8  =  (-  lkx^2^2  Jexp(-i/J2c) ,  (553) 

4.9  =  0  >  (554) 

4,io=0,  (555) 

4s,ii =0’  (556) 

4,12=0,  (557) 

4.1  =  (“  2^kxaz  )exp(i  a3c) ,  (558) 

4.2  =  (2/i3^a3)exp(-ia3c) ,  (559) 

4.3  =  (fhP?  ~  Vikx  )exp(i p3c) ,  (560) 

4.4  =  (M>2  ~  thkl  )exp(-i P3c)  ,  (561) 

4.5  =  i2^2kxa2  )exp(i  a2c) ,  (562) 

4.6  =  (-  2^2kxa2  )exp(-ia2c) ,  (563) 

4.7  =(-»  2P?  +  »2  kx  )exp(i£2c) ,  (564) 

4.8  =  (-  i“242  +  kx  )exp(-ij32c) ,  (565) 

4.9  =  0  >  (566) 

4.10  =  0  >  (567) 

4.11=0,  (568) 
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46,12-0.  (569) 

4,1=0,  (570) 

4,2  =  0,  (571) 

4,3=0,  (572) 

4,4=0,  (573) 

4.5  =(ia2)exp(ia2&)»  (574) 

4.6  =  (-ia2)expHa2&),  (575) 

4.7  =  6**)exp(i/M0»  (576) 

ArkWl-#).  (577) 

A79  =(-ia1)exp(ia1i),  (578) 

4,io  =(i«i)exp(-ia1&),  (579) 

4,n  =  (-  ikx  )exp(i fixb) ,  (580) 

4,12  =(-i^)exp(-i^&),  (581)  ' 

4,i=°.  (582) 

4,2=0,  (583) 

4,3=0,  (584) 

45,4=0*  .  (585) 

4s,5  =fe)exp(ia2^)*  (586) 

4.6  =fe)exp(-ia2fc),  (587) 

4.7  =  (-i/?2)exp(ij32&) ,  (588) 
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A, 8  =(i£2)exp(-ij32&). 

(589) 

4,9  =(-i^)exp(ia,fe), 

(590) 

4,io  =  (-i^)exp(-iai&). 

(591) 

4,n  =  OA  )exp(i/3,fc) , 

(592) 

4,12  =(-iA)exPHA^)> 

(593) 

4,i =  o » 

(594) 

O 

11 

<N 

(595) 

A>,3  ®  » 

(596) 

> 

II 

O 

(597) 

4,5  =(-«?/ 12  -2a? ^  -A2^)exp(ia2^) , 

(598) 

4,6  =  -2a?//2  -A^2  )exp(-ia2fc) , 

(599) 

4,7  =  (“  2^4iu2)exP(i^2&) . 

(600) 

4,8  =  faxPilh  )exp(-i/f2fc). 

(601) 

4,9  =  (a^  +  2a^i  +442)exp(ia^) , 

(602) 

4,io  =  (a?A,  +2a?^j  +  A,^)exp(-ia,&), 

(603) 

4,n  =(2MiA<i)exp(i0,&), 

(604) 

4i2  =  (-2^A^i)exp {-ifab). 

(605) 

> 

© 

II 

O 

(606) 

4o,2  _  0  > 

(607) 
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Ao,3  “O’ 


(608) 


Ao,4  =0,  (609) 

A  0,5  =  ("  2^2kxa2  )exp(ia2&) ,  (610) 

A  0,6  =  (2ft^«2)exPHa2i) ,  (611) 

A 0,7  =  (M?  ~  Vikl  )exp(i/?26) ,  (612) 

A 0,8  =  {^201  —  ^2kx  )exp(— i/^2^)  »  (613) 

A  0,9  =  (2Vikxa\  )exp(i«jfc) ,  (614) 

Ao,io  =(-2/ii^ai)exp(-ia,i>),  (615) 

Ao,n  =  (“  A*!#2  +  Vikx  )exp(i fab) ,  (616) 

Ao,i2  =  (~  +  Vikx  )exp(-ij0j&) ,  (617) 

Ai,i=°.  (618) 

4,2=0-  (619) 

Au=0.  (620) 

Ai,4=0,  (621) 

A  1,5  =°.  (622) 

Ai,6  =0.  (623) 

A  i,7=0.  (624) 

A  i,8  =  0  7  (625) 

A  1,9  =  (-  «'2A  -  2a, Vi  -  Aj  k\  )exp(i  ax  a) ,  (626) 


71 


A  uo  =  (-a?A  -  2a} ft  -  A^)exp(-ia,fl) , 

(627) 

A  ui  =  (-  2*,j8.0i  )exp(ijS,a) , 

(628) 

A  1,12  =  (2kxPiVi  )exp(-i/3,a) , 

(629) 

A  2,1  =  0  > 

(630) 

> 

k> 

i! 

o 

(631) 

A  2,3  =  o » 

(632) 

A  2, 4  =  0  , 

(633) 

A  2,5  =  0  » 

(634) 

A2,6  =  0  > 

(635) 

A  2,7  =  0  ’ 

(636) 

A2,8  =  0 , 

(637) 

A  2,9  =  (~  2Vikxa\  )exp(iaja) , 

(638) 

A  2,10  =  (2M\kxai  )exp(-iaja) , 

(639) 

A 2,11  =  ~  ftkx)exp(ifaa) , 

(640) 

A 2,1 2  =  (M*  -  Ai,^)exp(-i^a), 

(641) 

*1,1  =  MK’®)  * 

(642) 

*2.i  =B3(kx,(o), 

(643) 

*3,1  =  Cl(kx’®)  » 

(644) 

*4,i  =  B>-i{kx,(o) , 

(645) 
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x5i  =  A2(kx,co ) , 
x6  l  =  B2(kx,(o) , 

*7,1  =  ^2  ’ 

x81  =  D2(kx,co) , 

*9,1  =Al(kx,co), 

*10,1  =  Bi(kx,(o) , 

*11,1  =  > 

*12.1  =  A  >  ®) » 

&u  =  -2P[(kx,(0)exp(-iYd) , 
b2l  =  0 , 


b3l  =  0 , 

^4,1  =  ^ » 
^5,1  =  0 , 
^6,1  =  0  > 
&7 1  =  0 , 
\  l  =  °> 

&9  j  =  0 , 
^10,1  =  0  > 
^11,1  =  ^  ’ 


(646) 

(647) 

(648) 

(649) 

(650) 

(651) 

(652) 

(653) 

(654) 

(655) 

(656) 

(657) 

(658) 

(659) 

(660) 
(661) 
(662) 

(663) 

(664) 
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^12,1  -  0  • 


(665) 


Using  equations  (498)  -  (665),  the  solution  to  the  constants  A1 ,  Bx ,  C, ,  £>, ,  A2 ,  B2 ,  C2 ,  D2 , 
A3 ,  B3 ,  C3 ,  and  D3  can  be  found  by 


x  =  A_1b.  (666) 

Additionally,  the  transfer  function  between  the  normal  displacement  at  location  z  (when  z<b) 
and  the  structural  load  can  be  written  as 


T  (kx ,  co)  =  ■■  =  Al(kx,co)ia1  exp(i  a,  z)  - 

P/(kx,co)exp(-iyd) 

Bl(kx,co)  i  a,  exp(-i  a,  z)  + 

Cl(kx,  co)ikx  exp(i£,  z)  +  £>,  (kx ,  ft))i^  exp(-ift  z), 


or  (when  z>  b  and  z<c) 


(667) 


T (kx , £0)  =  ~7\  =A2(kx,co)itt2 exp(i a2z)~ 

Pl(kx,co)cxp(-iyd) 

B2(kx,co)ia2  exp(-i  a2  z)  + 

C2 {kx,co)ikx  exp(i p2z)  +  D2 ( kx , co)ikx  exp(-i )3,z), 


(668) 


of  (when  z>  c) 


T(k  ,co)  = 


Uz(z,kx,co) 


■  =  A3(kx,co)ia3  exp(ia3z)- 


P,(kx,co)exp(-iyd) 

B2(kx,co)  ia3  exp(-i  a3z)  + 

C3(kx,co)\kx  exp(i/?3z)  +  D3{kx,co)\kx  exp(-i/?3z). 


(669) 


The  transfer  function  between  the  normal  stress  at  location  z  (when  z<b)  and  the  structural  load 
can  be  written  as 
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T(-k*’(0)  p; (k~a)) exp(-iyd)  Al(k*,(0)(  2a'^  ^k2x)txp(ialZ)  ■ 

B1(kx,o))(-a^l-2affil  -A1^2)exp(-icc1z)  + 

Ci  (kx  >  ®)(-2^  )  exp(i£,  z)  +  Dx  (kx,o))(2kxplfil)exp(-iplz), 


(670) 


or  (when  z  >  b  and  z<c) 


T(k  ,(Q)  = 


T  (z,kx,(o) 


^2  (^x  >  ^)(  (^2^2  —0*2  f^2  ^2kx  "t” 


P/(£Jt,ty)exp(-iy<i) 

52 (fcx ,«)(- a2  A2  -  2a2V2  -  A2fc2  )exp(-ia,z)  + 

C2(kx,  co)(-2kx  P2p2)  exp(i  &  z)  +  D2  (kx ,  co)(2kx  P2p2)  exp(-i02  z),  (671) 


of  (when  z  >  c) 


T(kx,co)  = 


TL(z,ftx,ft>) 

P/(A:;t,(w)exp(-i7<i) 


=  A3(fc;c,£0)(-a32A3  -2a3Ju3  - A3fc 2 )exp(i a3z)  + 


B3(^:;c,<y)(-a3 A3  -2 a3p3  - A3£ 2 )exp(-i a3 z)  + 

Q  >  0))(~^kx  /l3p3 )  exp(i/?3  z)  +  D3(kx,  (0){2kx  /33p3 )  exp(-i/l3  z).  (672) 


Figure  23  is  a  plot  of  displacement  in  the  z-direction  (normal)  divided  by  normal  pressure 
versus  wavenumber  in  the  x-direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z-  d.  In  this 
example,  the  sum  of  the  thickness  of  the  three  plates  is  very  small  ( h  =  0.045  m).  The  solid  line 
represents  the  thick  plate  theory  using  three  plates  (equation  (669))  and  the  x’s  represent  the  thin 
plate  theory  (equation  (107))  using  h  =  0.045  m.  Additionally,  the  material  properties  of  the 
three  (thick)  plates  are  identical  so  that  a  direct  comparison  between  the  thick  and  thin  plate 
theories  can  be  made  using  only  a  single  thin  plate.  Parameters  used  to  formulate  this  model  are: 
bottom  of  plate  1  (a)  is  -0.045  m,  intersection  of  plates  2  and  1  ( b )  is  -0.030  m,  intersection  of 
plates  3  and  2  (c)  is  -0.015  m,  top  of  plate  3  ( d )  is  0  m,  Young’s  modulus  (E)  is  13.2e9  (l+0.03i) 

N/m2 ,  Poisson’s  ratio  (v)  is  0.30,  and  density  of  the  plate  (p)  is  1938  kg/m3 .  Only  positive 
values  of  wavenumber  are  shown,  as  the  function  is  symmetric  about  k  =  0.  Note  that  there  is 
almost  complete  agreement  between  the  thick  plate  theory  using  three  similar  thick  plates  and  the 
thin  plate  theory  using  one  plate  for  this  specific  example.  Figure  24  is  a  plot  of  the 
displacement  in  the  z-direction  (normal)  divided  by  normal  pressure  versus  wavenumber  in  the  x 
direction  at  a  forcing  frequency  of  500  Hz  evaluated  at  z  =  d.  In  this  example,  the  sum  of  the 
thickness  of  the  three  plates  is  very  small  (h  =  0.0032  m).  The  solid  line  represents  the  thick 
plate  theory  using  three  plates  (equation  (669))  and  the  x’s  represent  the  thin  plate  theory 
(equation  (107))  using  h  =  0.0030  m.  Additionally,  the  material  properties  of  the  three  (thick) 
plates  are  extremely  dissimilar  so  that  a  direct  comparison  between  the  thick  and  thin  plate 
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theories  can  be  made  using  only  a  single  thin  plate  with  the  thickness  of  the  stiff  thick  plate  (h  = 
0.0030  m).  Parameters  used  to  formulate  this  model  are:  bottom  of  plate  1  (a)  is  -0.0032  m, 
intersection  of  plates  2  and  1  ( b )  is  -0.0031  m,  intersection  of  plates  3  and  2  (c)  is  -0.0030  m,  top 

of  plate  3  00  is  0  m,  Young’s  modulus  of  plate  3  (£3)  is  13.2e9  (l+0.03i)  N/m2,  Young’s 
modulus  of  plates  2  and  1  (£2  and  £,)  is  13.2e4  N/m2,  Poisson’s  ratio  of  plates  1, 2,  and  3  (d) 

is  0.30,  and  density  of  plates  1, 2,  and  3  (p)  is  1938  kg/m3  .  Only  positive  values  of  wavenumber 
are  shown,  as  the  function  is  symmetric  about  k  =  0.  Note  that  there  is  almost  complete 
agreement  at  between  the  thick  plate  theory  using  three  dissimilar  thick  plates  and  the  thin  plate 
theory  using  one  plate  (matched  to  the  stiff  plate  parameters)  for  this  specific  example.  Figure 
25  is  a  plot  of  the  wavenumber-frequency  (kco)  surface  of  the  system  using  the  thick  plate 
equations  of  motion  (equation  (669))  evaluated  at  z  =  d.  The  first  plate  is  a  stiff  material  and  the 
parameters  of  this  plate  are:  bottom  of  the  plate  (a)  is  -0.08  m,  top  of  the  plate  (b)  is  -0.06  m, 

Young’s  modulus  of  the  plate  (£))  is  4.55el0  (l+0.03i)  N/m2,  Poisson’s  ratio  of  the  plate  (t>, ) 

is  0.30,  and  density  of  the  plate  ( pl )  is  7700  kg/m3 .  The  second  plate  is  a  soft  material  and  the 
parameters  of  this  plate  are:  bottom  of  the  plate  ( b )  is  -0.06  m,  top  of  the  plate  (c)  is  -0.01  m. 
Young’s  modulus  of  the  plate  (£2 )  is  le9(l+0.15i)  N/m2 ,  Poisson’s  ratio  of  the  plate  ( v2 )  is 

0.45,  and  density  of  the  plate  ( p2 )  is  1200  kg/m3 .  The  third  plate  is  a  stiff  material  and  the 
parameters  of  this  plate  are:  bottom  of  the  plate  (c)  is  -0.01  m,  top  of  the  plate  (d)  is  0  m, 

Young  s  modulus  of  the  plate  (E3)  is  13.2e9  (l+0.03i)  N/m2,  Poisson’s  ratio  of  the  plate  (t>3 )  is 

0.30,  and  density  of  the  plate  ( p3 )  is  1938  kg/m3 .  The  color  scale  to  the  right  of  the  plot  is  in 

decibels.  All  values  greater  than  -190  dB  are  displayed  as  -190  dB  and  all  values  less  than  -220 
dB  are  displayed  as  -220  dB. 

Figures  26,  27,  and  28  are  plots  of  stress  in  the  z-direction  (normal)  divided  by  normal 
pressure  versus  wavenumber  in  the  x-direction  at  forcing  frequencies  of  1000,  3000,  and  6000  Hz 
evaluated  at  z  —  ( b+c)/2 .  This  evaluation  point  is  the  middle  of  the  second  plate.  The  first  plate  is 
a  stiff  material  and  the  parameters  of  this  plate  are:  bottom  of  the  plate  (a)  is  -0.279  m,  top  of  the 

plate  ( b )  is  -0.203  m,  Young’s  modulus  of  the  plate  (£j)  is  2.07ell  N/m2,  Poisson’s  ratio  of  the 

plate  (t>, )  is  0.30,  and  density  of  the  plate  ( p, )  is  7830  kg/m3 .  The  second  plate  is  a  soft 
material  and  the  parameters  of  this  plate  are:  bottom  of  the  plate  ( b )  is  -0.203  m,  top  of  the  plate 
(c)  is  -0.076  m,  Young’s  modulus  of  the  plate  (£2)  is  le7(l+0.30i)  N/m2,  Poisson’s  ratio  of  the 

plate  (v2 )  is  0.40,  and  density  of  the  plate  ( p2 )  is  1000  kg/m3 .  The  third  plate  is  a  soft  material 
and  the  parameters  of  this  plate  are:  bottom  of  the  plate  (c)  is  -0.076  m,  top  of  the  plate  ( d)  is  0 
m.  Young’s  modulus  of  the  plate  (£3)  is  5e7  (1+0. 15i)  N/m2,  Poisson’s  ratio  of  the  plate  (u3 )  is 

0.45,  and  density  of  the  plate  ( p3)  is  1200  kg/m3 .  The  compressional  wavespeed  of  the  fluid 
(cf)  is  1500  m/s  and  the  density  of  the  fluid  (py)  is  1025  kg/m3 . 
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Phase  Angle  (degrees)  Magnitude  (dB) 


Figure  23.  Transfer  Function  of  a  Triple  Plate  With  Fluid  Load 
( Equation  (669))  Using  Similar  Material  Properties  at  500  Hz 


Phase  Angle  (degrees)  Magnitude  (dB) 


Figure  24.  Transfer  Function  of  a  Triple  Plate  with  Fluid  Load 
(Equation  (669))  Using  Dissimilar  Material  Properties  at  500  Hz 
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Figure  25.  Transfer  Surface  of  a  Triple  Plate  with  Fluid  Load 
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Figure  26.  Transfer  Function  of  a  Triple  Plate  with  Fluid  Load 
(Equation  (671))  at  1000  Hz 
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Phase  Angle  (degrees) 


Figure  27.  Transfer  Function  of  a  Triple  Plate  with  Fluid  Load 
( Equation  (671))  at  3000  Hz 
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Phase  Angle  (degrees) 


Figure  28.  Transfer  Function  of  a  Triple  Plate  with  Fluid  Load 
(Equation  (671))  at  6000  Hz 
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9.  SUMMARY 


The  transfer  functions  for  six  different  thick  plate  configurations  were  derived.  This 
included  a  single  thick  plate,  a  single  thick  plate  with  a  fluid  load,  a  double  thick  plate,  a  double 
thick  plate  with  a  fluid  load,  a  triple  thick  plate,  and  a  triple  thick  plate  with  a  fluid  load. 
Validation  comparisons  were  made  using  similar  transfer  functions  derived  with  the  thin  plate 
theory,  the  results  of  which  are  that  for  small  plate  thicknesses,  these  transfer  functions  are  nearly 
identical. 
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